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We extend the physics case for a new high-energy, ultra-high statistics neutrino scattering exper- 
iment, NuSOnG (Neutrino Scattering On Glass) to address a variety of issues including precision 
QCD measurements, extraction of structure functions, and the derived Parton Distribution Func- 
tions (PDFs). This experiment uses a Tevatron-based neutrino beam to obtain a sample of Deep 
Inelastic Scattering (DIS) events which is over two orders of magnitude larger than past samples. 
We outline an innovative method for fitting the structure functions using a parameterized energy 
shift which yields reduced systematic uncertainties. High statistics measurements, in combination 
with improved systematics, will enable NuSOnG to perform discerning tests of fundamental Stan- 
dard Model parameters as we search for deviations which may hint of "Beyond the Standard Model" 
physics. 
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A. NuSOnG: Precision Structure Functions and 
Incisive QCD Measurements 

The search for new physics at the "Terascale" - energy 
scales of ~ 1 TeV and beyond - is the highest priority for 
particle physics. 

NuSOnG is a proposed high energy, high statistics 
neutrino scattering experiment that can search for "new 
physics" from the keV through TeV energy scales via pre- 
cision electroweak and QCD measurements. 

During its five-year data acquisition period, the Nu- 
SOnG experiment could record almost one hundred thou- 
sand neutrino-electron elastic scatters, and hundreds of 
millions of Deep Inelastic Scattering (DIS) events, ex- 
ceeding the current world data sample by more than an 
order of magnitude. 

This experiment can address concerns related to ex- 
traction of structure functions and their derived Parton 
Distribution Functions (PDFs), investigate nuclear cor- 
rections, constrain isospin violation limits, and perform 
incisive measurement of heavy quarks. 
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II. DEEP INELASTIC SCATTERING AND 
PARTON DISTRIBUTION FUNCTIONS 



Experiment 


v DIS 

events 


v DIS 

events 


main 

target 


isoscalar 
correction 


CCFR 


0.95M 


0.17M 


iron 


5.67% [4J 


NuTeV 


0.86M 


0.24M 


iron 


5.74% [5] 


NuSOnG 


606M 


34M 


glass 


isoscalar 



Table I: Comparison of statistics and targets for parton distri- 
bution studies in NuSOnG compared to the two past highest 
statistics DIS neutrino scattering experiments. 



Obtaining a high quality model of the parton distri- 
bution functions in neutrino and antineutrino scattering 
is crucial to the NuSOnG electroweak measurements 0. 
NuSOnG will go a step beyond past experiments in ad- 
dressing the systematics of parton distribution functions 
(PDFs) by making high statistics measurements for neu- 
trino and antineutrino data separately. TableUshows the 
large improvement in statistics for NuSOnG compared 
to NuTeV and CCFR, the previous highest statistics ex- 
periments. Issues of uncertainties on the nuclear correc- 
tions are avoided by extracting PDFs on Si02 directly, in 
similar fashion to the NuTeV Paschos-Wolfenstein (PW) 
analysis. 

The differential cross sections for neutrino and antineu- 
trino CC DIS each depend on three structure functions: 
F2, xF 3 and Rl- They are given by: 
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where +(— ) is for ^(f) scattering. In this equation, x 
is the Bjorken scaling variable, y the inelasticity, and 
Q 2 the squared four-momentum transfer. The structure 
functions are directly related to the PDFs. 

The function xF 3 (x,Q 2 ) is unique to the DIS cross 
section for the weak interaction. It originates from the 
parity-violating term in the product of the leptonic and 
hadronic tensors. For an isoscalar target, in the quark- 
parton model, where s = s and c = c, 

xF 3 vN (x) = x (u(x) + d(x) + 2s(x) (2) 

—u{x) — d(x) — 2c(cc)) , 
xF» N (x) = xF£ N (x) - 4x (s(x) - c{x)) . (3) 

In past experiments, the average of xF 3 for neutrinos 
and antineutrinos has been measured. Defining xF 3 = 
\{xF% N + xF s " N ), at leading order in QCD, 

xF 3>LO = ^2 xqi(x,Q 2 ) - xq~(x,Q 2 ). (4) 

i—u.d.. 

To the level that the sea quark distributions have the 
same x dependence, and thus cancel, xF 3 can be thought 
of as probing the valence quark distributions. The differ- 
ence between the neutrino and antineutrino parity violat- 
ing structure functions, A(xF 3 ) = xF 3 N — xF 3 N , probes 
the strange and charm seas. (67/. Sec. IVII) 

The function F2(x,Q 2 ) appears in both the cross sec- 
tion for charged lepton (e or fx) DIS and the cross section 
for v DIS. At leading order, 

F 2 xo= e 2 {x qi {x,Q 2 )+xq-i(x,Q 2 )), (5) 

i— u,d.. 

where e, is the charge associated with the interaction. In 
the weak interaction, this charge is unity. For charged- 
lepton scattering mediated by a virtual photon, a is 
the fractional electromagnetic charge of the quark flavor. 
Thus Fg N and F% are analogous but not identical 
and comparison yields useful information about specific 
parton distribution flavors [3] and charge symmetry vio- 
lation as discussed below. In past neutrino experiments, 
F% and have been taken to be identical and an average 
F2 has been extracted, although this is not necessarily 
true in nuclear targets, as discussed below. 

Similarly, Rl(x, Q 2 ), the longitudinal to transverse vir- 
tual boson absorption cross-section ratio, appears in both 
the charged-lepton and neutrino scattering cross sections. 
To extract Rl from the cross section, one must bin in the 
variables x, Q 2 and y. This requires a very large data 
set. To date, the best measurements for Rl come from 



charged lepton scattering rather than neutrino scattering 
0. Therefore, neutrino experiments have used charged 
lepton fits to Rl as an input to the measurements of xF 3 
and F2 0- This, however, is just a matter of the statis- 
tics needed for a global fit to all of the unknown structure 
functions in x and Q 2 bins [§]. With the high statistics 
of NuSOnG, precise measurement of Rl will be possible 
from neutrino scattering for the first time. 

As an improvement on past experiments, the high 
statistics of NuSOnG allows measurement of up to six 
structure functions: F$ , F%, xF% , xF% , R L and R L . This 
is done by fitting the neutrino and antineutrino data sep- 
arately in X, y and Q 2 as described in Eq. ((J). The first 
steps toward fitting all six structure functions indepen- 
dently were made by the CCFR experiment [13], however 
statistics were such that only xF 3 u , xF 3 , and F 2 -average 
and i?-average could be measured, where the average is 
over v and v. A global fit of up to six structure functions 
in NuSOnG would allow separate parameterizations of 
the underlying PDFs which can account for the nuclear 
and isospin violation issues discussed below. 

In addition to fitting to the inclusive DIS sample, 
neutrino scattering can also probe parton distributions 
through exclusive samples. A unique and important case 
is the measurement of the strange sea through charged 
current (CC) opposite sign dimuon production. When 
the neutrino interacts with an s or d quark, it can produce 
a charm quark that fragments into a charmed hadron. 
The charmed hadron's semi-leptonic decay (with branch- 
ing ratio B c ~ 10%) produces a second muon of opposite 
sign from the first: 

+ N — ► n~ + c + X 

^ S + /!+ + Vfi . (6) 

Similarly, with antineutrinos, the interaction is with an 
s or d, 

Vf, + N — > /!+ + c + X 

n~ + 77 M . (7) 

The opposite sign of the two muons can be determined 
for those events where both muons reach the toroid spec- 
trometer. Study of these events as a function of the kine- 
matic variables allows extraction of the strange sea, the 
charm quark mass, the charmed particle branching ratio 
(B c ), and the Cabibbo-Kobayashi-Maskaka matrix ele- 
ment, \V cd \. 
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Figure 1: The assumed energy- weighted flux 

(Ed(j>/dE/10 6 POT) based on the NuTeV experiment 
in a) neutrino mode (left) and b) antineutrino mode (right), 
a) In neutrino mode the fluxes are ordered: upper (black), 
muon neutrino; middle (blue), electron neutrino and antineu- 
trino; lower (red), muon antineutrino. b) In antineutrino 
mode the fluxes are ordered: upper (red), muon antineutrino; 
middle (blue), electron neutrino and antineutrino; lower 
(black), muon neutrino. 



III. EXPERIMENTAL EXTRACTION OF 
STRUCTURE FUNCTIONS IN NUSONG 

A. Description of NuSOnG 

The NuSOnG detector was designed to be sensitive 
to a wide range of neutrino interactions from ^-electron 
scattering as well as ^-nucleon scattering, though this 
paper focuses mainly on the latter process. The design 
has been described in detail in Refs. 0, [H| and so we 
provide only a brief summary here. 

The neutrino beam would be produced via a high en- 
ergy external proton beam from a high intensity acceler- 
ator with energy of ~ 1 TeV. The Fermilab Tevatron is 
an existing example. The CERN SPS+ [H E3], which 
is presently under consideration because of its value to 
the LHC energy and luminosity upgrades and to a future 
beta beam, is a second example. 

The NuSOnG beam design will be based on the one 
used by the NuTeV experiment, which is the most recent 
high energy, high statistics neutrino experiment. The 
experiment would use 800 GeV protons on target fol- 
lowed by a quad-focused, sign-selected magnetic beam- 
line. The beam flux, shown in Fig. [H has very high neu- 
trino or antineutrino purity (~98%) and small v e con- 
tamination (~2%) from kaon and muon decay. Using an 
upgraded Tevatron beam extraction it is expected that 
NuSOnG could collect 5 x 10 19 protons/yr, an increase 
by a factor of 20 from NuTeV. With this high intensity, 
such a new facility would also produce a neutrino beam 
from the proton dump having a sizable fraction of tau 
neutrinos for study. 

The baseline detector design is composed of a 
fine-grained target calorimeter for electromagnetic and 
hadronic shower reconstruction followed by a toroid 



muon spectrometer to measure outgoing muon momenta. 
The target calorimeter will be composed of 2,500 2.5 cm 
x 5 m x 5 m glass planes interspersed with proportional 
tubes or scintillator planes. This gives a target which is 
made of isoscalar material with fine 1 /4 radiation length 
sampling. The detector will be composed of four target 
sections each followed by muon spectrometer sections and 
low mass decay regions to search for long-lived heavy neu- 
tral particles produced in the beam. The total length of 
the detector is ~200 m and the fiducial mass for the four 
target calorimeter modules will be 3 kiloton which is 6 
times larger than NuTeV or CHARM II. Figure [2] shows 
a simulated charged current event in the detector. 



B. Description of NuSOnG Calibration Beam 

The requirements for NuSOnG calibration beam would 
be similar to those of NuTeV. Tagged beams of hadrons, 
electrons, and muons over a wide energy range (5- 
200 GeV) would be required. The calibration beam will 
have the ability to be steered over the transverse face of 
the detector in order to map the magnetic field of each 
toroid with muons. Steering for hadrons and electrons 
would be less crucial than it was in NuTeV's case, but 
would still be useful. 

The calibration beam can be constructed with a sim- 
ilar design to NuTeV. Upstream elements were used to 
select hadrons, electrons, or muons. An enhanced beam 
of electrons was produced by introducing a thin lead radi- 
ator into the beam and detuning the portion of the beam 
downstream of the radiator. A radiator was also used 
in the nominal beam tune to remove electrons. Particle 
ID (a threshold Cerenkov and TRDs) was incorporated 
in the spectrometer and used to tag electrons when run- 
ning at low energy. A pure muon beam was produced by 
introducing a 7 m long beryllium filter in the beam as an 
absorber. 

The NuTeV calibration spectrometer determined in- 
coming particle momenta with a precision of better than 
0.3% absolute [3]. The NuSOnG goal for calibration- 
beam precision would be to measure energy scales to a 
precision of about 0.5%, and we demonstrate (in later 
text of this paper) that this can be improved with fits to 
neutrino data. 

For comparison, using the calibration beam, NuTeV 
achieved 0.43% precision on absolute hadronic energy 
scale and 0.7% on absolute muon energy scale (dominated 
by the ability to accurately determine the toroid map). 
Precise knowledge of the muon energy scale is especially 
important in order to achieve high measurement accuracy 
on the neutrino fluxes using the \ow-v method. For ex- 
ample, a 0.5% precision on muon energy scale translates 
into about a 1% precision on the flux. Both energy scales 
are important for precision structure function measure- 
ments, and were the largest contributions to structure 
function measurement uncertainties in NuTeV [l5| . 
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Process=91 , run=0, event=8 
Incoming: PDG=14, energy =446276 MeV 
vertex=(0.01 91 22,0.021 3528,404.54) cm 
# primaries=14. total particles=2342 
Em u=318897 MeV, theta=18.1429 mr, Etot-Emu=1 27379 MeV 
EmuFF=306724 MeV, theta=1 8.3256 mr 



Figure 2: A simulated muon neutrino, charged current event in the NuSOnG detector. 



Experimental Extraction of Structure Functions 
in NuSOnG 



The high statistics of the NuSOnG experiment makes 
it possible to extract the structure functions directly 
from the y-distributions within bins of (a;, Q 2 ). Previous 
lower-statistics high-energy neutrino experiments either 
extracted structure functions by comparing the number 
of v versus v events in an (x, Q 2 ) bin Q, or by extracting 
the cross-sections da/dy within the (x,Q 2 ) bin and fit- 
ting for the structure functions using Equation {!]) [TBI ] . 
Either method assumes a value for Rl = cfl/o't as mea- 
sured by other experiments flSl ]. and depends on a mea- 
surement of the strange sea from dimuon events [l7l.llgj]. 
With sufficient statistics, we can explore the possibility 
of measuring xF% (x, Q 2 ), xFg (x, Q 2 ) , F 2 (x,Q 2 ), and 
R (x, Q 2 ) from the same data llOj. 



Let us denote Eq. U) as a function of the struc- 
ture functions by da u ^ (xF 3 ,F 2 ,R), where the (x,Q 2 )- 
dependence is assumed and where the structure functions 
can be different for neutrinos and antineutrinos. A sam- 
ple of Monte Carlo events, N^c gen> ls generated using an 
assumed set of structure functions for the cross-section: 



ture functions in each 



n ). One can then fit for the struc- 
x,Q 2 ) bin by minimizing 
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events in an (x,Q 2 ,y) bin, is given by 



N MC,pred ( SF fit) = ^ d(J v ^ (xF^ F$ e " Rgen^ N MC,gen (S F gen) , (9) 
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Figure 3: Fractional change in number of events for two characteristic (x,Q 2 ) bins as a function of y. The fractional change 
comes from scaling the energy of each event by a factor of 1.005. 



^dati 1S tne number of v or v data events in the (x, y, Q 2 



bin, and N^\ gen 



is the number of Monte-Carlo 

■fit 



events generated in the (x,y,Q 2 ) bin. xF 3 lt (x, Q 2 ) , 
F2 (x, Q 2 ) , and R^ lt (x, Q 2 ) are the fit parameters in 
the ^-minimization of Eq. |(8|). In principle they can be 
fit separately for v and v structure functions. Here we 
will concentrate on the measurement ox f up to four sep- 
arate structure functions, AxF 3 (x,Q 2 ) = xF 3 (x, Q 2 ) — 
xFf(x,Q 2 ), xF™ 9 = (xF£ + zif)/2, F 2 (x,Q 2 ), 
and R (x, Q 2 ) where we assume that F2 (x, Q 2 ) and 
R (x, Q 2 ) are the same for neutrinos and antineutri- 
nos i.e. F 2 (x,Q 2 ) = F% (x,Q 2 ) = F% (x,Q 2 ) and 
R (x, Q 2 ) = R" (x, Q 2 ) = R~ (x, Q 2 ) . 



D. Fitting for AxF 3 

We have studied the extraction of the structure func- 
tion from the 600 million neutrino and 33 million anti- 
neutrino deep inelastic scattering events expected in the 
full NuSOnG data set. The dominant systematic error 
comes from the measurement of the muon momentum 
in the toroidal spectrometer. At NuTeV, the system- 
atic uncertainty was 0.7% and we assume NuSOnG will 
achieve 0.5%. Our studies are carried out by fitting the y- 
distribution in each x, Q 2 bin for F2, the average value of 
xF 3 = xF™ 9 = (xF£ + xFf)/2, AxF 3 = xF% - xFJ 7 ,and 
R. In the first set of studies, R(x, Q 2 ) is set equal to the 
measured value[16] and fits are done to the three struc- 
ture functions, F 2 , xF£ V9 , and AxF 3 . 

Our fitting procedure begins with a sample of Monte 
Carlo generated events, N gen {x, Q 2 , y), sampled from 
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Figure 4: Fractional uncertainty for the fit value of AxF$ in different x bins as a function of Q 2 . The fit is to multiple x and Q 2 
bins extracting the three structure functions, F2,xF£ V3 and AxFa. For each of the fits, a global set of energy scale parameters 
is also determined from the fit. The dotted lines show the fractional error for the NuTeV 2/^ measurement. 
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the CCFR structure functions and the nominal value for 
AxF 3 from NuTeV. We fit in bins of (x, Q 2 ) as a function 



of y and obtain the fit spectra by reweighting the original 
sample: 
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N° en (x,Q 2 ,y). 



r 



where the upper sign is for neutrinos and the lower for 
anti-neutrinos. In order to study the effects of the sys- 
tematic energy scale shift, we produce a Monte Carlo 
sample where the muon energy scale is shifted by 0.5%, 
E meas = i.005 J B* rue , for each event. The fractional 
change in the number of events in each bin due to the 
energy scale shift is shown in Fig. [3l 

This shifted event distribution, N shl f l {x, Q 2 , y), is 
then used to carry out a three parameter fit to Eq. [8] 
where F2, xF 3 vg , and AxF 3 are varied. Large shifts in 
AxF 3 result. For example, the shift from the input value 
in the (x, Q 2 ) = (0.08, \2.&GeV 2 ) bin is 19.01% and the 
shift in other bins is even larger. 

The effects of the energy scale uncertainty can be prac- 
tically eliminated by including energy scale shift param- 
eters in the fit. A muon energy scale change shifts the 
events in the various y-bins by an amount that is not con- 
sistent with that expected from changes in the structure 
functions. Therefore, fits to the y-distributions can iso- 
late the effects of an energy scale shift and significantly 
reduce the structure function uncertainty from this sys- 
tematic error. To estimate the systematic error reduc- 
tion for this technique, three additional energy scale pa- 
rameters are introduced in the fit to the y-distributions. 
These three parameters are used to produce an energy 
scale shift parameterization in each (x,Q 2 ,y) bin given 
by 



En 



Enscalel ~t~ F^ sca i e 2E^ 



F u,scalesF „ . 



The updated prediction for the number of events in a 
given (x,Q 2 ,y) bin is 

N pred( SF fit) = N MC,pred ( SF fit) 

+E liscale (N shift (x, Q 2 ,y) - N Ben (x, Q 2 ,y)) , 

(10) 

and the x 2 used in the minimization similar to Eq. [8] 
with the addition of pull terms associated with the three 
energy scale parameters 
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These pull terms correspond to an energy scale uncer- 
tainties of about 0.5% for muon energy values averaging 



between 50 and 70 GeV. This fitting technique renders 
the systematic error from the scale shift to be small in 
comparison with the statistical error. For example, in 
the bin (x, Q 2 ) = (0.275, 32 GeV 2 ) bin, the systematic 
error for AxF 3 is 0.3% while the statistical error is 10%; 
the value of the E usca i el parameter is also determined to 
about 10%. 



In the ultimate analysis, the fit will be carried out si- 
multaneously over all x and Q 2 bins with one set of en- 
ergy scale parameters. We have studied this using eight 
x bins and six to eight Q 2 bins. Figure [4] shows the frac- 
tional error on AxF 3 for different x bins as a function of 
Q 2 . In general, we believe NuSOnG can measure AxF 3 
over most of the {x,Q 2 ) range to better than 10%; in 
many cases around 3%. Typical values for NuTeV are 
shown in two x bins in Fig. [H Since more than one 
(x, Q 2 ) bin is being used to determine the energy scale 
shift parameters, the value of the -E^scaiei parameter can 
also determined to about 3% from these fits. 



Simulation studies have also been made to estimate the 
uncertainties associated with doing fits to extract the four 
structure functions, F2, xF^ V9 , Arri^and R. The proce- 
dure is the same as used for the three structure function 
fits where the \ 2 m Eq. Q1] is minimized simultaneously 
over a number of x and Q 2 bins with one set of energy 
scale parameters. In this case, the AxF 3 and Ri on g struc- 
ture functions can be determined to between 5% and 20% 
for most of the x and Q 2 range as shown in Figs. [5]and[6l 
The simulated Ri ong measurements are shown in Fig. [7| 
along with previous measurements. [l(| As indicated from 
this figure, the capabilities of the NuSOnG to measure 
Riong is much more precise that any previous experiment. 



In summary, due to the very high statistics of a Nu- 
SOnG type experiment, an almost complete set of struc- 
ture functions over a broad range of x and Q 2 can be 
extracted from the data without introducing theoreti- 
cal or experimental approximations. Further, systematic 
uncertainties that have limited the precision of previous 
structure function measurements can be eliminated by in- 
cluding fits to these uncertainties in the extraction proce- 
dure. We believe that with these techniques the structure 
function measurements will be statistics limited even for 
NuSOnG. 
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parameters is also determined. 
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IV. NUCLEAR EFFECTS 

Historically, neutrino experiments have played a major 
role in expanding our understanding of parton distribu- 
tion functions through hig h statistics experiments such as 
CCFR 01, NuTeV @,[IEl!1], and CHORUS [13]. How- 
ever, the high statistics extract a price since the large 
event samples require the use of nuclear targets - iron in 
the case of both CCFR and NuTeV and lead in the case 
of the Chorus experiment. The problem is that if one 
wants to extract information on nucleon PDFs, then the 
effects of the nuclear targets must first be removed. Nu- 
SOnG can provide key measurements which will improve 
these corrections. 

Charged lepton deep inelastic scattering has been mea- 
sured on a wide range of targets. The most simplistic ex- 
pectation for the structure functions might be that they 
would simply be given by an average of the appropriate 
number of proton and neutron results as in 



F 2 A (x,Q 2 ) = -rF 2 p (x,Q 2 ) + —^F?(x,Q 2 ). 



However, the results from a wide range of experiments 
show a much more complex behavior for the structure 
functions on nuclei. The typical behavior of the ratio 
of F^{x, Q 2 ) to F^x.Q 2 ) where d denotes a deuterium 
target shows four distinct regions as sketched in Fig. [HI 

At small x the ratio dips below one in what is called 
the shadowing region. At somewhat larger values of x 
the ratio rises above one in the antishadowing region. At 
still larger values of x the ratio again falls below one in 
the EMC region. Finally, as x approaches one, Fermi 
motion smearing causes a significant rise in the ratio. 

This behavior shows only a modest dependence on A 
for values above beryllium, with the shape remaining 
qualitatively the same and the amount of the suppression 
at x w 0.6 increasing slowly with \og(A). Furthermore, 
there is little, if any, observed dependence on Q 2 . These 
features are summarized nicely in the results shown in 
Ref. [H. 

The mechanisms of nuclear scattering have also been 
studied theoretically. These mechanisms appear to be 
different for small and large Bjorken x as viewed from the 
laboratory system. Bjorken x is defined as x = Q 2 /2Mv, 
where v and q are energy and momentum transfer to 
the target and Q 2 = q 2 — v 2 . The physical quantity 
which is responsible for the separation between large and 
small x regions is a characteristic scattering time, which 
is also known as Ioffe time (or length) rj = v/Q 2 
If tj is smaller than the average distance between bound 
nucleons in a nucleus then the process can be viewed as 
incoherent scattering off bound nucleons. This happens 
at larger x(> 0.2). 



A. Nuclear effects at small x 

We expect to find a difference between charged-lepton 
nucleus and neutrino nucleus scattering at small- a; be- 
cause the space-time pictures for the two processes are 
different in this region. The underlying physical mecha- 
nism in the laboratory reference frame can be sketched 
as a two-stage process. At the first stage, the virtual 
photon 7*, or W* or Z* in case of neutrino interactions, 
fluctuates into a quark-gluon (or hadronic) state. In the 
second stage, this hadronic state then interacts with the 
target. The uncertainty principle allows an estimate of 
the average lifetime of such hadronic fluctuation as 



m 2 + Q 2 xMm 2 + Q 2 ' y ' 

where m is invariant mass of hadrons into which the vir- 
tual boson convert, and M is the proton mass. The same 
scale t also determines characteristic longitudinal dis- 
tances involved in the process. At small x, r exceeds the 
average distance between bound nucleons. For this rea- 
son coherent multiple interactions of this hadronic fluc- 
tuation in a nucleus are important in the small-x kine- 
matical region. It is well known that the nuclear shad- 
owing effect for structure functions is a result of coherent 
nuclear interactions of hadronic fluctuations of virtual 
intermediate boson. 1 

For neutrino interactions which are mediated by the 
axial- vector current, the fluctuation time r is also given 
by Eq. [321 However, as was argued in Ref. [HI], the fluc- 
tuation and coherence lengths are not the same in this 
case. In particular, the coherence length is determined 
by the pion mass in Eq. [12] because of the dominance 
of off-diagonal transitions like a\N — > ttN in nuclear in- 
teractions. Since the pion mass is much smaller than typ- 
ical masses of intermediate hadronic states for the vec- 
tor current (ra p , m u , etc.), the coherence length L c of 
intermediate states of the axial current at low Q 2 will 
be much larger than L c of the vector current. A direct 
consequence of this observation is the early onset of nu- 
clear shadowing in neutrino scattering at low energy and 
and low Q 2 as compared with the shadowing in charged- 
lepton scattering. The basic reason for this earlier onset 
and different behavior in the transition region is the dif- 
ference in the correlation lengths of hadronic fluctuations 
of the vector and axial- vector currents. This is also illus- 
trated by observing that for a given Q 2 , the cross-section 
suppression due to shadowing occurs for much lower en- 
ergy transfer (v) in neutrino interactions than for charged 
leptons. 



1 For a recent review of nuclear shadowing see, e.g., (2^1 . 
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Figure 6: Fractional uncertainty for the fit value of R in different x bins as a function of Q 2 . The fit is to multiple x and Q 2 bins 
extracting the four structure functions, F2,xF^ V3 , AXF3, and R. For each of the fits, a global set of energy scale parameters is 
also determined. 
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B. Nuclear Effects in Neutrino Interactions 

As there has been no systematic experimental study 
of v and V nucleus interactions, one must then rely on 
theoretical models of the nuclear corrections. This is an 
unsatisfactory situation since one is essentially measuring 
quantities sensitive to the convolution of the the desired 
PDFs and unknown - or model dependent - nuclear cor- 
rections. 

As noted above, theoretically there are substantial dif- 
ferences between charged lepton and neutrino interac- 
tions on the same nucleus. There are other expected dif- 
ferences for neutrinos. For example, the relative nuclear 
shadowing effects for the structure function F3 is pre- 
dicted to be substantially different from that for F2 [25l |. 
This is because the structure function F3 describes the 
correlation between the vector and the axial- vector cur- 
rent in neutrino scattering. In terms of helicity cross 
sections, the structure function F3 is given by the cross 
section asymmetry between the left- and right-polarized 
states of a virtual W boson. It is known that such a dif- 
ference of cross sections is strongly affected by Glauber 
multiple scattering corrections in nuclei. [2a . |2~t1 . l28| 
This causes an enhanced nuclear shadowing effect for the 
structure function F3. 

It is important to experimentally address the question 
of nuclear effects in neutrino scattering so that the neu- 
trino data can be used in proton fits without bringing 
in substantial nuclear uncertainties. For example, in a 
recent analysis [2(| the impact of new neutrino data on 
global fits for PDFs was assessed. The conclusion reached 
in this analysis was that the uncertainties associated with 
nuclear corrections precluded using the neutrino data to 
constrain the nucleon PDFs. If NuSOnG can address 
these uncertainties, then the neutrino data can play an 
even more prominent role in the global fits to the proton 
PDF. 

Furthermore, nuclear effects are interesting in their 
own right. Parameterizations of nuclear PDFs on vari- 
ous targets exist in the literature. However, there is no 
universally accepted model which describes these nuclear 
corrections over the entire range of x from first principles. 
This makes it difficult to generalize the above behavior 
observed in charged lepton DIS to DIS with v or v beams. 
Models such as that in Ref. [30] exist, but to date there 
have been no high statistics studies of v or v DIS over a 
wide range of nuclear targets with which to test them. 

A study presented in Ref. examined the role of new 
lepton pair production data from E-866 and new neutrino 
DIS data from the NuTeV and CHORUS collaborations 
in global fits for nucleon PDFs. For the actual fitting of 
the PDFs it was necessary to include nuclear corrections 
for the neutrino and antineutrino cross sections and the 
model of Ref. (3(J was used. As a byproduct of that anal- 
ysis, it was possible to compare a reference fit, obtained 
without using data on nuclear targets, to the neutrino 
and antineutrino data in order to obtain an estimate of 
what the nuclear corrections should look like. This com- 



parison is shown in Fig. [9j 

This figure shows some results from Ref. [2§| in the 
form of "data/theory" averaged over Q 2 and presented 
versus x. The results are from a global fit but are plotted 
without the model-dependent nuclear corrections which 
were used in the fit (the neutrino data were not used in 
the reference fit.) It is notable that the overall pattern of 
deviations shown in Fig. [9] are, in general, similar to that 
seen in charged lepton DIS as sketched in Fig|8l However, 
the deviations from unity are perhaps smaller. At high 
x, the effect of Fermi smearing is clear. At moderate x 
the EMC effect is observable. It is interesting to note 
that there is no clear indication of a turnover at low x in 
the shadowing region for v data. Also, note the striking 
similarity between the v and V results. This appears to 
imply that the differences in the nuclear effects between 
neutrino and antineutrino DIS are small. As discussed 
later, when we consider AxF 3 and isospin violation, it is 
crucial to model differences in the nuclear effects between 
v and v scattering as a function of x. 

To make progress in understanding nuclear corrections 
in neutrino interactions, access to high statistics data on 
a variety of nuclear targets will be essential. This will 
allow the ^4-dependence to be studied as a function of 
both x and Q 2 , as has been done in charged lepton deep 
inelastic scattering. PDFs from global fits without the 
neutrino data can then be used to make predictions to be 
compared with the A-dependent v and V cross sections, 
thereby allowing the nuclear corrections to be mapped 
out for comparison with theoretical models. 

The primary target of NuSOnG will be Si02- However, 
we can investigate a range of A- values by replacing a few 
slabs of glass with alternative target materials: C, Al, 
Fe, and Pb. This range of nuclear targets would both 
extend the results of Minerva to the NuSOnG kinematic 
region, and provide a check (via the Fe target) against 
the NuTeV measurement. 

Given the NuSOnG neutrino flux, we anticipate 58fc v- 
induced and 30k v- induced CC DIS events per ton of ma- 
terial. A single ton would be sufficient to extract i<2(x) 
and xFz{x) averaged over all Q 2 ; a single 5 mx5 mx2.54 
cm slab of any of the above materials will weigh more 
than that. The use of additional slabs would permit fur- 
ther extraction of the structure functions into separate 
(x, Q 2 ) bins as was done in the NuTeV analysis, at the 
potential expense of complicating the shower energy res- 
olution in the sub-detectors containing the alternative 
targets; this issue will be studied via simulation. 

Table HI] shows that two 50-module stacks would be 
sufficient to accumulate enough statistics on alternative 
nuclear targets for a full structure-function extraction for 
each material. However, for basic cross-section ratios in 
x, a single slab of each would suffice. 
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Figure 7: Extracted values of R (labeled NuSOnG) from the four structure function fits as compared to previous 
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Measuring Nuclear Effects with the Minerva 
and NuSOnG Detectors 



The Minerva experiment will also be studying neutrino 
induced nuclear effects and will be starting its initial 
physics run in early 2010. To study nuclear effects in 
Minerva, a cryogenic vessel containing liquid helium (0.2 
ton fiducial mass)will be installed upstream of the Min- 
erva detector. Within the Minerva detector, solid car- 
bon, iron and lead targets will be installed upstream of 
the pure scintillator active detector. The total mass is 
0.7 ton of Fe, 0.85 ton of Pb, 0.4 ton of He and some- 
what over 0.15 ton of C. Since the pure scintillator active 
detector essentially acts as an additional 3-5 ton carbon 
target (CH), the pure graphite (C) target is mainly to 



Material 


Mass of 
2.54 cm slab (tons) 


Number of slabs needed 
for NuTeV-equivalent statistics 


C 


1.6 


33 


Al 


1.9 


27 


Fe 


5.5 


10 


Pb 


7.9 


7 



Table II: Alternative target materials for cross-section analy- 
sis 



check for consistency. For a run consisting of 4.0 x 10 20 
POT in the NuMI Low Energy (LE) beam and 12 x 10 20 
POT in the NuMI Medium Energy (ME) beam, Minerva 
would collect over 2 M events on Fe, 2.5 M events on Pb, 
600 K on helium and 430 K events on C as well as 9.0 M 
events on the scintillator within the fiducial volume. 

Studying nuclear effects with the NuSOnG detector 
will involve fewer nuclear targets but considerably more 
statistics on each. In addition, the much higher energy 
of the incoming neutrinos with NuSOnG means a much 
wider kinematic range of study. In particular, NuSOnG 
will have a much higher Q 2 for a given low- a; to study 
shadowing by neutrinos and will be able to measure the 
shadowing region down to much smaller x for the same 
Q 2 range as Minerva. A significant addition to the study 
of nuclear effects with neutrinos would be the addition of 
a large, perhaps active ("Bubble Chamber"), cryogenic 
target containing hydrogen or deuterium. With the in- 
tense NuSOnG neutrino beam, a significant sample of 
neutrino-hydrogen and neutrino-deuterium events could 
provide the normalization we need to further unfold nu- 
clear effects in neutrino-nucleus interactions. 



V. QCD FITS 

The extraction of up to six structure functions from 
the cross sections of neutrino and anti-neutrino DIS dis- 
cussed so far (c/, Eq. ((TJ)) has been completely model- 
independent relying only on some fundamental principles 
such as Lorentz-invariance of the cross section and gauge- 
invariance of the hadronic tensor which is expanded in 
terms of the structure functions which parameterize the 
unknown hadronic physics. 

More can be said about the structure functions in 
QCD. While it is still not possible to accurately com- 
pute the x-dependence of the structure functions from 
first principles, QCD allows us to derive renormalization 
group equations (RGEs) which relate the structure func- 
tions at different (perturbative) scales Q. Note that the 
structure functions at the scale Q can be directly related 
to structure functions at a different scale Qo (see, e.g., 
Eqs. (5.58) and (5.76) in [HI]). However, it is more conve- 
nient to work in the QCD-improved parton model where 
the RGEs governing the scale-dependence of the parton 
distribution functions (PDFs) are the familiar DGLAP 
evolution equations; these can also be used to compute 
the structure functions at Q given the PDFs at that scale. 
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[32, |33j, |34| Furthermore, this approach has the crucial 
advantage that the universal PDFs allow us to make pre- 
dictions for other observables as well. In addition to 
the Q-dependence, the QCD calculations provide certain 
(approximate) relations between different structure func- 
tions as will be visible from the parton model expressions 
below. 

In this section we will discuss the analysis of the cross 
section data within the framework of the QCD-improved 
parton model. Already in the past, high statistics mea- 
surements of neutrino deeply inelastic scattering (DIS) on 
heavy nuclear targets (NuTeV, ...) have attracted much 
interest in the literature since they provide valuable in- 
formation for global fits of PDFs (HEEl- 

Due to the weak nature of neutrino interactions, the 
use of nuclear targets is unavoidable; this complicates 
the extraction of free nucleon PDFs, because model- 
dependent corrections must be applied to the data (cf 
Sec. lIVp . Of course, these same data are also useful 
for extracting the nuclear parton distribution functions 
(NPDFs) and for such an analysis no nuclear correc- 
tion factors are required. Conversely, the NPDFs can 
be utilized to compute the required nuclear correction 
factors within the QCD parton model [32} • Similar to 
proton PDFs, universal nuclear PDFs are needed for the 
description of many processes with nuclei in the initial 
state. This involves physics at other neutrino experi- 
ments, heavy ion colliders (RHIC, LHC), and a possible 
future electron-ion collider (EIC). 

The NuSOnG experiment will have two orders of mag- 
nitude higher statistics than the NuTeV and CCFR ex- 
periments (over an extended kinematic range), and so 
it will be possible to study small effects such as the 
strangeness asymmetry with better precision, or to estab- 
lish for the first time isospin violation in the light quark 
sector. Better understanding these effects is relevant for 
improving the extraction of the weak mixing angle in a 
Paschos-Wolfenstein type analysis. 



A. PDFs 

NuSOnG will perform measurements on different nu- 
clear targets. The PDFs for a nucleus (A, Z) are con- 
structed as 



A 



In the following discussion we take into account devia- 
tions from isospin symmetry, a non- vanishing strangeness 
asymmetry and the possibility to have non-isoscalar tar- 
gets. For this purpose we introduce the following linear 
combinations of strange quark PDFs: 



s+> A = s A + s A , s-> A = s A -s A 



(14) 



where the strangeness asymmetry is described by a non- 
vanishing PDF s~ . Note however that we continue to 



assume s p / A = s n l A and s p / A = s n / A . Also, we neglect 
any possible charm asymmetry, i.e., we use c A = c A such 
that c~' A — c A - c A — and c + - A = c A + c A = 2c A . 

Deviations from isospin symmetry can be parameter- 
ized in the following way: 

8v?J A = u p J A - d n J A , 8d p J A = d p J A - u n J A , (15) 

Su p/A = u p/A - d n/A , Sd p/A = d p/A - u n/A . (16) 

These definitions allow us to write the PDFs in a way 
which makes deviations from isoscalarity and isospin 
symmetry manifest: 



2u 



- [vPj A + dPj A - 5d p J A ] - 
A[u p J A - d p J A + 5d p / A ] , 
2d A = [u%/ A + <Pj A -5ul' A ] + 
A[u p J A - d p / A - Su p / A ] , 



(17) 



(18) 



2u A = [vP' A + d p > A - 5d p ' A ] - 

A[u p / A - d p ' A + Sd p/A ] , (19) 

2d A = [vP' A + dP/ A -8v?/ A ] + 

A[u p / A - d p l A - 8u p / A ] , (20) 

where A = (N — Z)/A parameterizes the deviation from 
isoscalarity. We have written Eqs. (fl7| - l[20|) so that the 
RHS is expressed explicitly in terms of proton PDFs 
and the four <5-terms {5u? M , 5d p J A , 8u p / A , 8d p l A }; the 
(5-terms vanish individually if isospin symmetry is pre- 
served. 

B. Structure functions 

The structure functions for a nuclear target (A, Z) are 
given by 

F A (x,Q) = j F[/ A ( X ,Q) + F:' A {x,Q) (21) 

such that they can be computed in next-to-leading order 
as convolutions of the nuclear PDFs with the conven- 
tional Wilson coefficients, i.e., generically 



F A (x,Q)=Y / C*k®f A . 



(22) 



In order to discuss which information can be extracted 
from a high statistics measurement of neutrino and anti- 
neutrino DIS cross sections we briefly review the parton 
model expressions for the 6 structure functions. For sim- 
plicity, we first restrict ourselves to leading order, neglect 
heavy quark mass effects (as well as the associated pro- 
duction thresholds), and assume a diagonal CKM matrix. 
In our numerical results, these effects are taken into ac- 
count. 

The neutrino-nucleus structure functions are given by 
(suppressing the dependence on x and Q 2 ): 

F" A = d A + s A + u A + c A + ... , (23) 

' jA - 2xF£ A , (24) 

1 - c A + . . .] . (25) 



F» A 



2 [d A + s A 
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Figure 10: Structure function ratio Rl for neutrino and anti- 
neutrino nucleon (TV — (p + n)/2)) scattering at Q 2 = 20 
GeV 2 . The solid and dashed lines show NLO results obtained 
with the GRV98 PDFs [H], while the dotted line shows the 
LO result of Eq. (T3TT) . 



The structure functions for anti-neutrino scattering are 
obtained by exchanging the quark and anti-quark PDFs 
in the corresponding neutrino structure functions: 



The differential cross section in Eq. (Q]) can be written 



as: 



with K 



da 
dxdy 



K[A + B(l-y) 2 + Cy 2 } 



(32) 



at, me 



2^(i+q»/M^)2 . A = F * ± xF 3, B = F 2 T xF 3 , 
and C = 2x q¥ Fi — Fl where the upper sign refers to 
neutrino and the lower one to anti-neutrino scattering. 
This form of da shows that the (anti-)neutrino cross sec- 
tion data naturally encodes information on the four struc- 
ture function combinations F% A ±xF% A and F£ A ± xF% A 
in separate regions of the phase space. In addition, at 
large y the structure function combination C contributes. 
However, to good accuracy C v = C u so that C drops out 
in the difference of neutrino and anti-neutrino cross sec- 
tions. 



Assuming s J 



and 



c , the structure func- 



tions F% A and F£ A constrain the valence distributions 



d A = d A 
sea S A := 



d A , ug = u 
u A + d A + s 

l -F» A 

x 

X -Fl A 



Furthermore, we have 



u 

7-A 



A and the flavor-symmetric 
via the relations: 



2 [d A + ^ 
2 [u A + X J 



(33) 
(34) 



Explicitly this gives 



Ft 



-F» A {q^q}. (26) 



X 

1 



pA 



F 



jpuA 
^3 



4{d A + s A ) 



4(<f 



* ' s A ) 



(35) 
(36) 



I 

F 



PA 



u A + c A + d* 



PA 



PA 



2xF{ 

2 k 



PA 



-d 1 



(27) 
(28) 
(29) 



The longitudinal structure function can be obtained with 
the help of the following relation: 



Since we constrain the strange distribution utilizing the 
dimuon data, the latter two structure functions are useful 
to separately extract the d A and d A distributions. 

For an isoscalar nucleus we encounter further simplifi- 
cations. In this case, u A = d A and u A = d A =: q A which 
implies u A = d A = 
distributions are {u", q' 



r 2 F 2 " 



Hence, the independent quark 
i , . . .}. In par- 



s A , c A 



. 2jy,p ticular, we have F% = F% for an isoscalar target such 

F% A , (30) that our original set of 6 independent structure functions 
^ reduces to 3 independent functions (say F% A , F£ A , F% A ) 



where r 2 = 1 + 4x 2 M 2 /Q 2 . Finally, it is customary to 
introduce the ratio of longitudinal to transverse structure 
functions: 



Fl A 

?pA 



r 2 F uA 

Ha 



Ax 2 M 2 



2xF» A 2xF{ 



(31) 



under the approximations made. 

In a more refined analysis, allowing for a non- vanishing 
strangeness asymmetry and isospin violation we can eval- 
uate the non-singlet structure function AF^ = F% A — 
F% A with the help of the relations in Eqs. fT7 |) - IpO ]) : 



AF 2 " 



Similar equations hold for anti-neutrino scattering. As 
can be seen, in leading order R V L = R V L . As is shown 
in Fig. Hoi also in NLO, the differences between R V L and 
R V L are tiny such that the difference between these two 
functions can be neglected in the following discussion. 



2xs-' A + x 5dl' A -x Su p J A 

+Ax[2uP/ A - 2dPj A + 5dPj A - 5vPj A ] . (37) 



For a nuclear isoscalar target (Z 
this expression simplifies to 



N 



AFi 



2xs 



-.A 



Sdf/ A 



A/2, A = 0) 
x 5vPj A . (38) 
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As one can see, Ai^ will be small and sensitive to the 
strangeness asymmetry and isospin violating terms for 
the valence quarks. 

The difference of the neutrino and anti-neutrino cross- 
sections provides, in principle, access to this quantity: 



2„vA 



d 2 <T 



•%„T>A 



d 2 a 



dxdy dxdy 



K[AF%+xFZ 

+{l-yf{AF^-xF^)} (39) 



F£ A . 



with F% = F% A 

It should be noted, however, that in a global fit to 
extract structure functions we do not make direct use 
of these equations [the (1 — y) 2 -dependence] but simply 
perform a x 2_ana ly s i s of all neutrino and anti-neutrino 
cross section data. 



VI. ISOSPIN (CHARGE SYMMETRY) 
VIOLATION AND AxF 3 

The question of isospin violation is central to the 
PW electroweak measurement. In the NuTeV analysis, 
isospin symmetry was assumed. As discussed in Ref. [3], 
various models which admit isospin violation can pull the 
NuTeV sin#jy measurement toward the Standard Model. 
However it would take significantly larger isospin viola- 
tion to bring NuTeV into agreement with the rest of the 
world's data. Better constraints of isospin violation will 
be crucial to the interpretation of the NuSOnG results. 

When we relate DIS measurements from heavy targets 
such as SjgFe (used in NuTeV) or §§ 7 Pb (Chorus) back 
to a proton or isoscalar target, we generally make use 
of isospin symmetry where we assume that the proton 
and neutron PDFs can be related via a u *-> d inter- 
change. While isospin symmetry is elegant and well mo- 
tivated, the validity of this exact charge symmetry must 
ultimately be established by experimental measurement. 
There have been a number of studies investigating isospin 
symmetry violation [H, S3, S3, S3, S3, SI EH] ; therefore, 
it is important to be aware of the magnitude of potential 
violations of isospin symmetry and the consequences on 
the extracted PDF components. For example, the naive 
parton model relations are modified if we have a violation 
of exact p <-> n isospin-symmetry, or charge symmetry vi- 
olation (CSV); e.g., u n {x) dP{x) and u p (x) ^ d n (x). 

It is noteworthy that a violation of isospin symmetry is 
automatically generated once QED effects are taken into 
46, 



account m, 143, |48 . This is because the photon couples 
to the up quark distribution u p (x) differently than to the 
down quark distribution d n (x). These terms can be as 
much as a few percent in the medium x range, see e.g. 
Fig. 1 in Ref. [H]. 

Combinations of structure functions can be particu- 
larly sensitive to isospin violations, and NuSOnG is well 
suited to measure some of these observables. For exam- 
ple, residual u, ^-contributions to AxF 3 = xFg — xFg 
from charge symmetry violation would be amplified due 



to enhanced valence components {u v (x), d v (x)}, and be- 
cause the d u transitions are not subject to slow- 
rescaling corrections which suppress the s — » c contri- 
bution to AxF 3 0. Here the ability of NuSOnG to 
separately measure xF% and xF 3 over a broad kinematic 
range will provide powerful constraints on the sensitive 
structure function combination AxF 3 . 

Separately, the measurement of AF2 = 
^Ff c (x,Q 2 ) - Fi fc (x,Q 2 ) in Charged Current 
(CC) W ± exchange and Neutral Current (NC) 7/Z 
exchange processes can also constrain CSV [43|; because 
NuSOnG will measure F^ on a variety of targets, this 
will reduce the systematics associated with the heavy 
nuclear target corrections thus providing an additional 
avenue to study CSV. 

In the following, we provide a detailed analysis of CSV 
which also investigates the various experimental system- 
atics associated with each measurement. We shall find it 
is important to consider all the systematics which impact 
the various experimental measurements to assess the dis- 
criminating power. 



A. AxF-j, and Isospin Violations 

We recall the leading-order relations of the neutrino 
structure function F3 on a general nuclear target: 



1 



*3 (*) 



d A + s A 



u A -c A 



-F* A (x) = u A + c A -d A -s A 



(40) 
(41) 



where A represents the nuclear target A — {p,n,d, ...}, 
and the represent higher-order contributions and 

terms from the third generation {&, t} quarks. Note that 
to illustrate the general features of these processes, we 
use a schematic notation as in Eq. l(40|) and Eq. (|4T|); 
for the numerical calculations, the full NLO expressions 
are employed including mass thresholds, "slow-rescaling" 
variables, target mass corrections, and CKM elements 
where appropriate. 

For a nuclear target A we can construct AxF A as: 



AxF( 



xFt 



2xA 



A - 2xc+' A 



p/ A _ (P>/ A 
+ 0(a s ) 



A 



-5I J 



(42) 



where O (as) represents the higher order QCD correc- 
tions, and the isospin violations are given by 51 : 



SI A = Sd— Su + 6d- Si 



(43) 



For a flux- weighted linear combination of F£ and F£, 
terms proportional to the stran g e q uark asymmetry can 
enter Eq. d42j), cf. Refs. [ISSaSil]- For a sign-selected 
v jv beam as for NuTeV or NuSOnG, this complication is 
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not necessary. We have defined s ±,A (x) = [s A (x)±s A (x)] 
and c^ix) = [c A {x)±c A {x)\. 

In the limit of isospin symmetry, all four terms on 
the RHS of Eq. lj43"|) vanish individually. For a nuclear 
isoscalar target, Z = N = A/2, we can construct AxF^ 
from the above: 



with the definitions: 



AxF, 



A xFg A 



2xs 



+- A -2xc + ' A +x5I A +0{a s ) 
(44) 

Note in Eq. f42|) that for a nuclear target A which is 
close to isoscalar we have Z ~ N such that the up and 
down quark terms are suppressed; this is a benefit of 
the NuSOnG glass (S1O2) target which is very nearly 
isoscalar. More specifically, for SiC>2 we have Z(0) = 8, 
Z(Si) = 14, m(O) = 15.994, m(Si) = 28.0855. Using A = 
Z + N we have (N-Z)/A = (A-2Z)/A for the prefactor 

in Eq. (02J| which yields (N - Z)/A 0.000375 for O 

and (N -Z)/A~ 0.00304 for Si. 

In Eq. @g) the PDFs {u p / A ,d p / A , ...} represent quark 
distributions bound in a nucleus A. With a single nuclear 
target, we can determine the CSV term 6I A for this spe- 
cific A; measurements on different nuclear targets would 
be required in order to obtain the A dependence of 5I A 
if we need to scale to a proton or isoscalar target. 

Thus, an extraction of any isospin violation 51 A re- 
quires a careful separation of these contributions from 
the strange, charm, and higher order terms. Theoretical 
NLO calculations for AxF$ are available; thus the O (as) 
corrections can be addressed. Additionally, NuSOnG can 
use the dimuon process (uN — » n + [i~X) to constrain the 
strange sea. 

In conclusion we find that while this is a challenging 
measurement, NuSOnG's high statistics measurement of 
AxF^ should provide a window on CSV which is rela- 
tively free of large experimental systematics. We empha- 
size that AXF3 may be extracted from a single target, 
thereby avoiding the complications of introducing nu- 
clear corrections associated with different targets. This 
is in contrast to the other measurements discussed below. 
However, if we desire to rescale the SI A effects to a differ- 
ent nucleus A, then multiple targets would be required. 



F, 



CC,A 



= 5 [JT 



pNC,A _ plA 

In Eq. (|45j) . the first term is proportional to (N — Z)/A 
which vanishes for an isoscalar target. The second and 
third terms are proportional to the heavy quark distri- 
butions s + ' A and c + ' A . The next term is the CSV con- 
tribution which is proportional to SI given in Eq. (|43|) . 
It is curious that this has the same form as the CSV con- 
tribution for AxF% of Eq. ([4"2"]l. Finally, the last term 
represents the higher-order QCD corrections. 

While the character of the terms on the LHS of Eq. |44|) 
and Eq. (|45|) are quite similar, the systematics of mea- 
suring AF 2 may differ substantially from that of AxF 3 . 
For example, the measurement of AF2 requires the sub- 
traction of structure functions from two entirely differ- 
ent experiments. The CC neutrino-nucleon data are ex- 
tracted from heavy nuclear targets (to accumulate suffi- 
cient statistics); as such, these data are generally subject 
to large nuclear corrections so that the heavy targets can 
be related to the isoscalar N — \(p + n) limit. Con- 
versely, the NC charged-lepton-nucleon process proceeds 
via the electromagnetic interaction. Therefore sufficient 
statistics can be obtained for light targets including H 
and D and no large heavy target corrections are neces- 
sary. Therefore, we must use the appropriate nuclear 
correction factors when we combine F^ c and F^ c , and 
this will introduce a systematic uncertainty. 

Separately, the heavy quark production mechanism is 
different in the CC and NC processes. Specifically, in the 
CC case we encounter the process s + W + — > c where the 
charm mass threshold kinematics must be implemented. 
On the other hand, the NC process is c + 7 — > c which is 
proportional to the charm sea distribution and has dif- 
ferent threshold behavior than the CC process. Even 
though the charm production process is modeled at NLO, 
the theoretical uncertainties which this introduces can 
dominate precision measurements. 



B. Measurement of AF 2 = A Ff c (x, Q 2 ) - Ff c (x, ( 



Other Measurements of CSV 



A separate determination of CSV can be achieved using 
the measurement of F 2 in CC and NC processes via the 
relation: 



AF 2 



5 pec 
18 2 

1 (N — Z) 
6*— 4— 



(uP/A _ d P/A^ + ^ P /A _ d p,A^ 



X s 



\x)-\xc + < A {x) + \x^6I A 



+ 0(a s ) 



(45) 



We very briefly survey other measurements of CSV in 
comparison to the above. 

The measurement of the lepton charge asymmetry in 
W decays from the Tevatron can constrain the up and 
down quark distributions [i^, [H3|- In this case, the ex- 
traction of CSV constraints is subtle; while isospin sym- 
metry is not needed to relate p and p, this symmetry is 
typically used in a global fit of the PDFs to reduce data 
on heavy targets to p. 

In the limit that all the data in the analysis were from 
proton targets, CSV would not enter; hence this limit 
only arises indirectly from the mix of targets which enter 
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a global fit. At present, while much of the data does 
come from proton targets (HI, ZEUS, CDF, DO), there 
are some data sets from both p and d (BCDMS, NMC, 
E866), and some that use heavier targets (E-605, NuTeV) 
[29L l5l| . Thus, an outstanding question is if CSV were 
present, to what extent would this be "absorbed" into a 
global fit. The ideal procedure would be to parameterize 
the CSV and include this in a global analysis. While this 
step has yet to be implemented, there is a recent effort 
to include the nuclear corrections as a dynamic part of a 
global fit [al. 

Additionally, NMC measures F% / F$ data which has an 
uncertainty of order a few percent [52]]. There are also 
fixed- targ et Drell-Yan experiments such as NA51 [45[ and 
E866 [53j] which are sensitive to the ratio d/u in the range 
0.04 < x < 0.27. We will soon have LHC data (pp) to add 
to our collection, thus providing additional constraints in 
a new kinematic region. 



D. Conclusions on Charge Symmetry Violation 

NuSOnG will be able to provide high statistics DIS 
measurements across a wide x range. Because the target 
material (SiCh) is nearly isoscalar, this will essentially 
allow a direct extraction of the isoscalar observables. 

AxF 3 is one of the cleaner measurements of CSV in 
terms of associated experimental systematic uncertain- 
ties as this measurement can be extracted from a single 
target. The challenge here will be to maximize the event 
samples. 

The measurement of AF 2 is more complicated as this 
must combine measurements from both CC and NC ex- 
periments which introduces nuclear correction factors 
[STI . Hij]. Since NuSOnG will provide high statistics 
^'"measurements for a variety of A targets, this will 
yield an alternate handle on the CSV and also improve 
our understanding of the associated nuclear corrections. 

The combination of these measurements, together with 
external constraints, will yield important information on 
this fundamental symmetry. 



VII. MEASUREMENTS OF THE HEAVY 
QUARKS 

A. Measurement of the Strange Sea 

Charged current neutrino-induced charm production, 
(v/v)N — > /i + /i~A, proceeds primarily through the sub- 
processes W + s — > c and W~s — > c (respectively), so 
this provides a unique mechanism to directly probe the 
s(x) and s(x) distributions. Approximately 10% of the 
time the charmed particles decay into jj, + X, adding a 
second oppositely signed muon to the CC event's final 
state. These "dimuon" events are easily distinguishable, 
and make up approximately 1% of the total CC event 
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Figure 11: NuTeV measurement of xs~ (x) vs x at Q 2 = 
16GeV 2 . Outer band is combined errors, inner band is with- 
out B c uncertainty. 

sample. Hence, the recent hig h statistics dimuon mea- 
surements [17L fTsl . [55l . I56L [53] play an essential role in 
constraining the strange and anti-strange components of 
the proton. On NuSOnG, the dimuon data will be used 
in the same manner. 

Distinguishing the difference between the s(x) and s(x) 
distributions, 

xs~(x) = xs(x) — xs(x), (46) 

is necessary for the PW style analysis. This analysis is 
sensitive to the integrated strange sea asymmetry, 

S~ = [ s~(x)dx, (47) 
Jo 

through its effect on the denominator of the PW ratio, as 
has been recognized in numerous references [H, [H, 0, 

EH El). 

The highest precision study of s~ to date is from 
the NuTeV experiment (H, [H. The sign selected 
beam allowed measurement of the strange and anti- 
strange seas independently, recording 5163 neutrino- 
induced dimuons, and 1380 antineutrino induced dimuon 
events in its iron target. Figure [TT] shows the fit for asym- 
metry between the strange and anti-strange seas in the 
NuTeV data. 

The integrated strange sea asymmetry from NuTeV 
has a positive central value: 0.00196 ± 0.00046 (stat) 
±0.00045 (syst) iaooio? (external). In NuSOnG, as in 
NuTeV, the statistical error will be dominated by the an- 
tineutrino data set and is expected to be about 0.0002. 
The systematic error is dominated by the ir and K decay- 
in-flight subtraction. This can be addressed in NuSOnG 
through test-beam measurements which will allow a more 
accurate modeling of this background, as well as applying 
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the techniques of CCFR to constrain this rate |G j.l6a.l6a|. 
We expect to be able to reduce this error to about 0.0002. 
The combination of these reduces the total error by about 
10%, because the main contributon comes from the ex- 
ternal inputs. 

The external error on the measurement is dominated 
by the error on the average charm semi-muonic branching 
ratio, B c : 



(48) 



where <fi ls the neutrino flux in energy bins, fa is the 
energy dependent production fraction for each hadron, 
and B^-i is the semi-muonic branching ratio for each 
hadron. In the NuTeV analysis, this is an external input, 
with an error of about 10%. To make further progress, 
this error must be reduced. 

Fig. [32] shows the world measurements of B c , taken 
from references [13, \M, H3, H3, HI, H3, Measuring 
B c directly requires the capability to resolve the indi- 
vidual charmed particles created in the interaction. The 
best direct measurements are from emulsion. This kind 
of measurement has been performed in past experiments 
(E531, Chorus) using emulsion detectors (69L l7(il|. where 
the decay of the charmed meson is well tagged. Since 
the cross section for charmed meson production is en- 
ergy dependent, it is important to make a measurement 
near the energy range of interest. The NuTeV strange sea 
asymmetry study used a re-analysis of 125 charm events 
measured by the FNAL E531 experiment [69] in the en- 
ergy range of the NuTeV analysis (E v > 20 GeV) . B c has 
also been constrained through indirect measurement via 
fits. 

For NuSOnG, our goal is to reduce the error on B c 
using an in situ measurement on glass by at least a factor 
of 1.5. One method is to include B c as a fit parameter 
in the analysis of the dimuon data. The unprecedentedly 
high statistics will allow a fit as a function of neutrino 
energy for the first time. Dimuons from high x neutrino 
DIS almost exclusively result from scattering off valence 
quarks, such that the dimuon cross section in that region 
isolates B c from the strange sea. In dimuon fits, the 
assumption is then taken that B c _ v = B c _„, B c may be 
measured directly from the dimuon data. 

Unfortunately, antineutrino charm production is not 
well measured by past experiments. This leads to con- 
cerns about the assumption that B c - V = B c _„. An 
example of a potential source of difference in neutrino 
and antineutrino mode, consider that vn — * /i~A c has no 
analogous reaction in the antineutrino channel. 

These arguments provide the motivation for including 
a high resolution target/tracker in the NuSOnG design 
that can directly measure the semileptonic branching ra- 
tio to charm in both v and v running modes. There 
are two feasible detector technologies. The first is to use 
emulsion, as in past experiments. This is proven technol- 
ogy and scanning could be done at the facility in Nagoya, 
Japan. The second is to use the NOMAD-STAR detector 
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Fig ure 12: World measurements of B c . See refs. [171. [5a. [57,. 



[71|,[72| or a similar detector. This is a 45 kg silicon vertex 
detector which ran in front of the NOMAD experiment. 
The target was boron carbide interleaved with the silicon. 
This detector successfully measured 45 charm events in 
that beam, identifying D + , D° and D s . A similar detec- 
tor of this size in the NuSOnG beam would yield about 
900 v events and 300 v events. This has the advantage 
of being a low-Z material which is isoscalar and close in 
mass to the Si02 of the detector. 



B. Strange Quark Contribution to the Proton Spin 

An investigation of the strange quark contribution to 
the elastic vector and axial form factors of the proton 
is possible in NuSOnG, by observing NC elastic and 
CC quasi-elastic scattering events; namely vp — > vp and 
vn — > u~p events in neutrino mode, and vp — > vp and 
vp — > u + n events in antineutrino mode. The motivation 
for making this measurement comes from a number of 
recent (and not so recent) studies in proton structure. 

Over the last 15 years a tremendous effort has been 
made at MIT-Bates, Jefferson Lab, and Mainz to mea- 
sure the strange quark contribution to the vector form 
factors (that is, the electromagnetic form factors) of the 
proton via parity-violating electron scatte ring fro m p ro- 
tons^deuterons, and 4 He [300S0lHlSIS£llsl 
Is3,ls3|- The technique is to observe the parity- violating 
beam spin asymmetry in elastic scattering of longitudi- 
nally polarized electrons from these unpolarized targets; 
this asymmetry is caused by an interference between the 
one-photon and one-Z exchange amplitudes [84] . As a 
result, the weak neutral current analog of the electro- 
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magnetic form factors of the proton may be measured 
and this gives access to the strange quark contribution. 
This worldwide experimental program will soon be com- 
plete. The results available to date (from global analy- 
ses [85l.[86l.[87|) indicate a small (and nearly zero) contri- 
bution of the strange quarks to the elastic electric form 
factor, G S E ; this is not surprising, as the total electric 
charge in the proton due to strange quarks is zero. At 
the same time, these same data point to a small but likely 
positive contribution of the strange quarks to the elastic 
magnetic form factor, G S M , indicating a small positive 
contribution of the strange quarks to the proton mag- 
netic moment. Due to the prominent role played by the 
Z-exchange amplitude, these experiments are also sensi- 
tive to the strange quark contribution to the elastic axial 
form factor, which is related to the proton spin structure. 

It is now well established by leptonic deep inelastic 
scattering experiments that the spins of the valence and 
sea quarks in the proton together contribute about 30% of 
the total proton intrinsic angular momentum oiTi/2. The 
strange quark contribution is estimated to be about -10% 
in inclusive DIS (an analysis which makes use of SU(3)- 
flavor symmetry) [13], but is found to be approximately 
zero in semi-inclusive DIS (an alternative analysis which 
makes no use of SU(3) but needs fragmentation functions 
instead) [!§]. A recent global analysis [§3] which made 
use of both inclusive and semi-inclusive DIS and which 
allowed for the possibility of SU(3)-flavor violation found 
no need in the data for any violation of SU(3) and in- 
dicated a small negative contribution of strange quarks 
to the proton spin. In the deep inelastic context, the 
contribution strange quarks make to the proton spin is 
encapsulated in the helicity-difference strange quark par- 
ton distribution function, 

As(x) = s~*(x) — s* - (x) 

where s~*(x) [s <_ (x)] is the probability density for find- 
ing a strange quark of momentum fraction x with its spin 
parallel [anti-parallel] to the proton spin. The axial cur- 
rent relates the first moment of this parton distribution 
function to the value of the strange quark contribution 
to the elastic axial form factor of the proton [9l|, G 8 a , at 
Q 2 -0: 

/ dxAs(x) = G S A {Q 2 =0). 
Jo 

The strange quark contribution to the elastic axial form 
factor can be measured by combining data from neutrino 
NC elastic scattering from the proton with data from 
parity- violating elastic ep scattering [92] • In this way 
the strange quark spin contribution to the proton spin 
can be measured in a completely independent way using 
low-Q 2 elastic scattering instead of high-Q 2 deep inelas- 
tic scattering. An analysis done using this method [85| 
indicates that G S A may in fact be negative at Q 2 = but 
this conclusion is not definitive due to the limitations of 
the currently available neutrino data. 



Since the neutrino experiments will undoubtedly be 
carried out on nuclear targets (perhaps carbon or argon) , 
then the extraction of the properties of the proton from 
these data needs to be done with care. Recent theoretical 
investigations point to the idea of measuring the ratio of 
NC to CC yields; nuclear effects appear to largely cancel 
in this ratio [§3| , leaving behind the ratio that would have 
been obtained on nucleon targets. 

The only available data on neutrino NC elastic scat- 
tering is from the BNL E734 experiment [94] ; the uncer- 
tainties reported from that experiment are considerable 
and limit the preciseness of any extraction of G S A based 
on them. If NuSOnG can provide more precise measure- 
ments of NC elastic scattering extended to lower Q 2 then 
the promise of this analysis technique can be fulfilled. 

C. Measurements of the Charm Sea: 

1. Charm Production 

We can also study the charm sea component of the 
proton which can arise from the gluon splitting pro- 
cess g — ► cc producing charm constituents inside the 
proton. [9|| [IB, [s3] In a measurement complementary to 
the above strange sea extraction, the charm sea, c(x,/x), 
can be measured using the following process: 

f)x + C — > Vfj,+ c 

S + /Lt + + 

In this process, we excite a constituent charm quark in 
the proton via the NC exchange of a Z boson; the fi- 
nal state charm quark then decays semi-leptonically into 
sp^v^. We refer to this process as Wrong Sign Muon 
(WSM) production as the observed muon is typically the 
opposite sign from the expected v^d — > fi~u DIS process. 
For antineutrino beams, there is a complementary pro- 
cess Dp + c —* i/p + c with a subsequent c — > s + + 
decay with yields a WSM with respect to the conven- 
tional v^u — ► fi + d process. Here, the ability of NuSOnG 
to have sign-selected beams is crucial to this measure- 
ment as it allows us to distinguish the secondary muons, 
and thus extract the charm-sea component. 

In the conventional implementation of the heavy quark 
PDFs, the charm quark becomes an active parton in the 
proton when the scale [i is greater than the charm mass 
m c \ i.e. f c (x,fj,) is nonzero for [i > m c . Additionally, we 
must "rescale" the Bjorken x variable as we have a mas- 
sive charm in the final state. The original rescaling proce- 
dure is to make the substitution x — > x(l+m 2 /Q 2 ) which 
provides a kinematic penalty for producing the heavy 
charm quark in the final state. [H] As the charm is pair- 
produced by the g —* cc process, there are actually two 
charm quarks in the final state — one which is observed in 
the semi-leptonic decay, and one which decays hadroni- 
cally and is part of the hadronic shower. Thus, the ap- 
propriate rescaling is not x — > x(l + m 2 /Q 2 ) but instead 
x — > x — - T (l + 4m 2 /Q 2 ); this rescaling is implemented 
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in the ACOT-% scheme, for example. [H Eol The 
factor (1 + 4m 2 . /Q 2 ) represents a kinematic suppression 
factor which will suppress the charm process relative to 
the lighter quarks. 

The differential cross section for NC neutrino scatter- 
ing is 



da 
d£,dy 



{up 



G pMjyEv 



R 2 z(Q 2 



M 



E„ 



£c(£,/x), 



where ql = £3 — Qf. sin 9w, 9r = —Qi sin 9w, and for 
charm t 3 = 1/2 and Q c = 2/3. The factor R Z {Q 2 ) = 
1/(1 + Q 2 /Mz) arises from the Z-boson propagator. The 
corresponding result for the anti-charm is given with the 
substitutions <7l «-> gu and c <-> c. 

In the limit we can neglect the M-^fE v term we 
have the approximate expressions for the total cross 
section: 1 96] 
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Figure 13: Integrated momentum fractions J Q x fi{x,Q) of 
charm (upper curve) and bottom (lower curve) PDFs (in per- 
cent) vs. Q in GeV. Both the quark and antiquark contri- 
butions are included. Horizontal lines at 0.5% and 1.0% are 
indicated as this is the typical size of postulated intrinsic con- 
tributions. 



and 



a(up 



uc) 



a{up — > uc) 



G 2 F M^E U 



G 2 F M]\iE v 



(0.129) C 



(0.063) C 



(49) 



(50) 



with C = X 



l min £ c(£, M ) d£ and C = £ c(£, /x) df. We 
take C = a:(l + 4m^/Q 2 ) and Cm in = m 2 J{2M N u). 

We will be searching for the WSM signal compared to 
the conventional charged-current DIS process; therefore 
it is useful to benchmark the rate for WSM production 
by comparing this to the the usual charged-current DIS 
process, 



massive charm cross section is down a factor of ~ 0.005 
compared to the total inclusive cross section. As the 
muon from the NC charm process is a secondary muon, 
we must additionally fold in the semi-leptonic branching 
ratio B c ~ 10%, and the acceptance factor of observ- 
ing the secondary muon in the detector (A,~20%).[l7j 
Combining the relevant factors, we estimate the rate for 
NC charm production is approximately a factor of 10~ 4 
as compared to the CC DIS process. Thus, for an antici- 
pated design of 600M u^ CC events, one would expect on 
the order of 60K NC charm events. This estimate is also 
consistent with a direct scaling from the NuTeV result of 
Ref. [H. 



2. Backgrounds 



da 
dx dy 



I - v \ G f MnE u 2 2 

{up -> p x) = R W \Q ) x 

X [q{x) + (1 - y) 2 q{x)} (51) 



with Rw{Q 2 ) = 1/{1 + Q 2 /Mw)- We can again integrate 
over x and y to obtain an estimate of the total cross 
section in terms of the integrated PDFs as in Eq. lj49j) 
and Eq. ([50]) : 



a{uN 



li-X) 
1 



G f MnE v 



RwiQ 2 ) x 



U + D + 2S+-{U + D + 2C) 



(52) 



where {U,D,S} are defined analogously to C, and we 
have used N = |(p + n) for an isoscalar target. 

The relative rate for NC charm production is deter- 
mined by the above factors together with a ratio of inte- 
grated PDFs. For a mean neutrino energy of 100 GeV, the 



Extrapolating from investigations by CCFR [96j, and 
NuTeV [97|, the dominant background for the measure- 
ment of the charm sea comes from u^ contamination. 
In these studies, it was determined that by demanding 
E v i s > 100 GeV, the background rate could be reduced 
to 2.3 x 10 -4 . Other background processes include u e in- 
duced dilepton production, mis-identified dimuon events, 
and NC interactions with a n/ K decay in the hadron 
shower; these processes contribute approximately an ad- 
ditional 1.5 x 10~ 4 to the background rate. As compared 
to CCFR and NuTeV, the NuSOnG design has a number 
of improvements such as lower mass density for improved 
shower measurement; hence, comparable background re- 
ductions should be achievable. 



3. Intrinsic Charm 

In the above discussion we have assumed that the 
charm component of the proton arises perturbatively 
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from gluons splitting into charm quark pairs, g — > cc; 
in this scenario the charm PDF typically vanishes at 
scales below the charm mass (f c (x, fi < m c ) = 0), and for 
fi > m c all the charm partons arise from gluon splitting. 

There is an alternative picture where the charm quarks 
are taken to be intrinsic to the proton; in this case there 
are intrinsic charm partons present at scales \i < in c . 
For \i > m C) the charm PDF is then a combination of 
this "intrinsic" PDF and the "extrinsic" PDF component 
arising from the g — ► cc process. 

A number of analyses have searched for an intrinsic 
charm component of the proton, and this intrinsic com- 
ponent is typically constrained to have an int egrated m o- 
mentum fraction less than a percent or two [l02l Il03| . 

In Figure[l3]we display the integrated momentum frac- 
tion, J Q x fi(x, /it), for charm and bottom as a function of 
H due to the "extrinsic" PDF component arising from 
the g — > cc or g — > bb process. These momentum frac- 
tions start from zero at the corresponding quark mass, 
and increase slowly as the partonic components pick up 
momentum from the gluon splitting process. 

If we are searching for an additional intrinsic compo- 
nent with a momentum fraction of ~ 1%, we will be 
most sensitive to such a component in the threshold re- 
gion where the "intrinsic" component is not overwhelmed 
by the "extrinsic" contribution. In this regard, NuSOnG 
is well suited to search for these intrinsic terms as it will 
provide good statistics in the threshold region. Measur- 
ing the charm production process described above, Nu- 
SOnG can attempt to extract the charm PDF as a func- 
tion of the /Li scale, and then evolve back to /i = m c . 
Three outcomes are possible: 

1. f c (x, fx = m c ) < 0, which would imply the data are 
inconsistent with the normal QCD evolution. 2 

2. f c (x,/i = m c ) = 0, which would imply the data is 
consistent with no intrinsic charm PDF. 

3. / C (a;,/Li = m c ) > 0, which would imply the data is 
inconsistent with an intrinsic charm PDF. 

By making accurate measurements of charm induced pro- 
cesses in the threshold region, NuSOnG can provide a 
discriminating test to determine which of the above pos- 
sibilities is favored. Hence, the high statistics of NuSOnG 
in the threshold region are well suited to further constrain 
the question of an intrinsic charm component. 

VIII. SUMMARY AND CONCLUSIONS 

The NuSOnG experiment can search for "new physics" 
from the keV through TeV energy scales. This article 



has focused mainly on the QCD physics which can be 
accessed with this new high energy, high statistics neu- 
trino scattering experiment. During its five-year data ac- 
quisition period, the NuSOnG experiment could record 
almost one hundred thousand neutrino-electron elastic 
scatters and hundreds of millions of deep inelastic scat- 
tering events, exceeding the current world data sample 
by more than an order of magnitude. 

With this wealth of data, NuSOnG can address a wide 
variety of topics including the following. 

• NuSOnG can increase the statistics of the Elas- 
tic Scattering (ES) and Deeply Inelastic Scattering 
(DIS) data sets by nearly two orders of magnitude. 

• The unprecedented statistics of NuSOnG allow the 
possibility to perform separate extractions of the 
structure functions: {F% , xF% , R V L , F 2 P , xF% ,B? L }. 
This allows us to test many of the symmetries 
and assumptions which were employed in previous 
structure function determinations. 

• NuSOnG will help us to disentangle the nuclear 
effects which are present in the PDFs. Furthermore, 
this may help us address the long-standing tensions 
between the NC charged- lepton and CC neutrino 
DIS measurements. 

• High precision NuSOnG measurements are sensi- 
tive to Charge Symmetry Violation (CSV) and 
other "new physics" processes. Such effects can sig- 
nificantly influence precision Standard Model pa- 
rameter extractions such as sinQy/ ■ In particular, 
AxF 3 is a sensitive probe of both the heavy quark 
components, and CSV effects. 

• NuSOnG dimuon production provides an excep- 
tional probe of the strange quark PDFs, and the 
sign-selected beam can separately study s{x) and 
s(x). Additionally, NuSOnG can probe the s-quark 
contribution to the proton spin. 

• The high statistics of NuSOnG may allow the mea- 
surement of the charm sea and an method to prove 
the intrinsic-charm content of the proton. While 
this is a difficult measurement, the NuSOnG kine- 
matics allow the measurement of charm-induced 
processes in the threshold region where the "intrin- 
sic" character can most easily be discerned. 

While the above list presents a very compelling physics 
case for NuSOnG, this is only a subset of the full range 
of investigations that can be addressed with this facility. 
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